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This work studies the mathematics of selected techniques for coloring fractal 
images. The classic escape time algorithm is extended by adding the concepts of 
coloring, palette and index functions. The coloring function is evaluated for each 
pixel of an image, whereas palette and index functions map this value to the final 
RGB color. In addition to good performance, also smoothness and the possibility 
to adjust the visual appearance are desirable characteristics of coloring functions. 

The Smooth Iteration Count coloring is the smooth equivalent of the classic 
Iteration Count coloring. Its continuity and smoothness are studied thoroughly. 
The Smooth Iteration Count coloring exhibits small discontinuities at the iter- 
ation boundaries of a fractal, and a method for calculating an upper bound for 
the magnitude of the discontinuity is presented. 

Average colorings are a family of coloring functions that use the decimal part 
of the smooth iteration count to interpolate between average sums. In addition 
to linear interpolation, a smooth Catmull-Rom spline interpolation method can 
be used. 

The Triangle Inequality Average and Curvature Average colorings are pre- 
sented and analyzed as examples of branching average colorings. Both colorings 
exhibit a similar tree-like branching structure, and the three mathematical prop- 
erties that lead to this kind of structure are presented. A new coloring called the 
Stripe Average is introduced based on the behavior of the Triangle Inequality 


Average coloring. 
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Chapter 1 
Introduction 


This work studies the mathematics of selected algorithms that are used to color 
fractal images. These algorithms were developed to be instruments for fractal 
art. Javier Barrallo and Damien Jones describe the emergence of fractal art and 


the role of coloring algorithms in the following way [2]. 


During the 1980s, fractal enthusiasts began exploring fractals for their 
artistic merit, not for their mathematical significance. While mathe- 
matics was the tool, the focus was art. As the fractal equation itself 
was the most obvious mathematical element, fractal artists experi- 
mented with new equations, introducing hundreds of different fractal 
types. By carefully choosing parameters to refine form, color, and 


location, these explorers introduced the concept of fractal art. 


After 1995, few new major fractal types have been introduced. This 
is because the newest innovations in fractal art do not come from 
changing the fractal equation, but from new ways of coloring the 
results of those equations. As these coloring algorithms move from 
simple to complex, fractal artists are often returning to the simpler, 
classical fractal equations. With the increased flexibility these so- 
phisticated algorithms provide, there is even more room for personal 


artistic expression. 


Countless coloring algorithms have been developed. The development pro- 
cess of a coloring algorithm usually combines experimentation to mathematical 


understanding and analysis of the visual results. However, thorough mathe- 
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matical understanding is not required in order to develop interesting and useful 
algorithms. 

This work aims to gain better understanding of the selected coloring algo- 
rithms, or colorings. Their visual appearance is explained based on mathemati- 
cal analysis. This analysis also led to the development of a new coloring, namely 
the Stripe Average coloring to be introduced in Section 4.5. 

Although this work aims to present its contents rigorously and accurately, 
some of the results are expressed in a somewhat vague form. This is inherent to 
the applied nature of the subject. To name an example, the words ”small” and 
* distant” are used loosely without an exact definition. 

The reader is assumed to be familiar with calculus and complex analysis. 
Fractals are explained in the extent necessary to introduce the main concepts 
and notation, whereas details in the theory of chaotic dynamical systems from 
which fractals emerge fall outside the scope of this work. 

The contents are organized as follows: Fundamental concepts and definitions 
necessary for the later chapters are introduced in Chapter 2. The classic escape 
time algorithm for calculating fractal images is extended to suit the needs of 
coloring functions. Heuristic criteria to evaluate the quality of a coloring are also 
presented. 

Chapter 3 presents two continuous colorings that extend the discontinuous 
Iteration Count coloring. A detailed discussion of the Smooth Iteration Count 
coloring and a proof of its smoothness are presented. It exhibits small disconti- 
nuities at the iteration boundaries of a fractal, and a method for calculating an 
upper bound for the magnitude of the discontinuity is presented. 

Chapter 4 first introduces average colorings followed by a discussion of linear 
and smooth Catmull-Rom spline interpolation methods. The Triangle Inequality 
Average and Curvature Average colorings are described and analyzed as examples 
of average colorings. Based on the behavior of the Triangle Inequality Average 
coloring, the Stripe Average coloring is then introduced. 

Appendix A illustrates the construction process of a fractal image using the 
colorings discussed in the earlier chapters. 

The thesis was carried out under the supervision of Professor Goran Hognas. 


Kerry Mitchell provided some helpful comments. 


Chapter 2 
Fundamental Concepts 


This chapter first introduces fundamental definitions that are the foundation for 
the following chapters. The calculation of images of divergent fractals is then 
discussed. 

For convenience, the presentation of definitions is adapted for the needs of 
this work. For example, definitions are presented as confined to the complex 


plane instead of an arbitrary metric space. 


2.1 Iterations and Orbits 


The notation and definitions in this section are adapted from [3]. 
The n-fold composition of f with itself is denoted by f(z). That is, 


where f is applied to z n times. The points f(z), f?(z), f?(z), ... are called the 
iterates of z. Repeatedly applying the function to the previous iterate is called 
iteration. Given an initial value zo and a function f, the nth iterate f"(zo) is 


denoted z,. 


Definition 2.1 (Orbit). Given a function f : C+ C, the (forward') orbit of a 
point z € C is the set 


O(z) = {2, F(2), F(2), ++}. 


'Tt is not necessary to separate forward and backward orbits in this context because only 


the forward orbit is of interest for this work. 
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The orbit O(z) is said to be periodic if, for a given f and z, there exists a 
number n € N such that f"(z) = z. An orbit contains infinitely many points if 
it is not periodic. An orbit O(z) escapes if the iterates converge to infinity. 

When calculating fractals on a computer, a finite representation of an orbit is 
needed. This is done by limiting both the magnitude and the maximum number 
of elements to be contained in an orbit. Define the constants M and Nyax as the 
bailout value and the maximum number of iterations, respectively. Furthermore, 


define the set Cy,,,,. as 
Nmax +1 


Cree Cr 
k=1 
Definition 2.2 (Truncated Orbit). Let an orbit 


OB =I OF i xa} 
and constants M and Ninax € N be given. Let N denote the smallest nonnegative 
integer for which |f%(z)| > M, and define N = min{N, Nmax}. The truncated 
orbit Or(z) € Cy,,,, is the set 
Or(z) = 12 F (2); F(Z), Caer] icone 


In addition to the initial point z, the truncated orbit contains N of the first 
Nmax iterates z,. The number of elements in a truncated orbit is N+ 1. It is 
noteworthy that N < Nya, and N is dependent of both the function f and the 
initial value z. 

The following inequality describes the magnitude of the two last iterates zy_1 


and zy of a truncated orbit with respect to the bailout value M. 


Proposition 2.1 (Bailout Inequality). Assume 1 < N < Ninax and consider the 
points zn-1 € Or(z) and zy € O7(z). The bailout inequality is the inequality 


|zn-1| <M< |zn]. (2.1) 


According to Definition 2.2, N is the smallest integer for which M < |zy}. 
Thus |zy_1| < M and the Bailout Inequality follows. 


2.2 Fractals 


As an example, consider the function f : Cr C, 


fig)=2 +e (2.2) 
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where the constant p € N, p > 2, and the seed c€ C. The function f(z) defines 
the dynamical system 

Be Se AE, (2.3) 
This system is known to exhibit chaotic behavior; depending on the value of the 
constant c and the initial value zg, the iterates behave differently. The intricate 
nature of this function can be illustrated by plotting some characteristics of the 


orbits when varying c or zp over the complex plane. 


2.2.1 Coloring and Palette functions 


An image is represented as an m xX n matrix of discrete points referred to as 
pixels. To each pixel in the matrix is associated an RGB color [4, p. 572]. Let zo 
denote the point corresponding to the position of a pixel in the complex plane. 
To calculate the RGB color of a pixel in a fractal image, the truncated orbit 


Or(2o) is first calculated. Then a coloring function is evaluated. 


Definition 2.3 (Coloring Function). A coloring function is a function 


u:Cn,,.. 0 R (2.4) 


max 


that maps a truncated orbit to the real line. 


The number of elements N +1 in the argument of u is variable. Implementa- 
tions of coloring functions usually include a loop that is executed for all elements 
of Or(zo) (cf. Section 2.2.2). The coloring u(Or(zo)) is often studied as a func- 
tion of the point zg, whereas certain coloring functions only use the last iterate 
zy of Or(z0). This is denoted u(zy). 

A color index function I(u) maps the value u of the coloring function to a 
color index. The color index is either contained in or clamped to the interval 
[0,1]. The color palette function P(I) is then used to map the color index to the 
RGB color space. The red, green and blue components of an RGB color can be 


represented as a point in (0, 1]°. 


Definition 2.4 (Palette Function). A palette function P : R ++ R® maps the 
color index I to the RGB color space. Its domain is {0,1] and range [0,1]°. The 


domain can be extended to the real line by using the decimal part of I. 


The color index function is merely used as a means to adjust the visual 


appearance of a coloring without the need to change the palette function. 
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2.2.2 Calculating Images of Fractals 


Iterating the function defined by equation 2.3, the algorithm presented below 
renders an image of a Julia set [1, p. 258]. This algorithm extends the escape 
time algorithm presented in [1, p. 252] by introducing coloring, index and palette 
functions. 

The following operations are carried out for each pixel in the m x n matrix 


representing the image. 


1. Set z to correspond to the position of the pixel in the complex plane. 


2. Calculate the truncated orbit by iterating the formula z,, = f(z,_1) starting 


from zp until either 


e |z,| > M, or 


en= Nese, 
where Nymax is the maximum number of iterations. 


3. Using the coloring and color index functions, map the resulting truncated 


orbit to a color index value. 


4. Determine an RGB color of the pixel by using the palette function. 


If |z,| > M, the pixel is an outside point. Otherwise, if m = Nmax, the point is 
inside. 

The Mandelbrot set can be calculated by setting z9 = 0 and c to correspond to 
the position of the pixel in the complex plane in step 1. In addition to the system 
2.2 mainly used in this work, many other functions exhibit chaotic behavior. 
Also various other types of fractals exist and are calculated with modified or 


completely different algorithms, but are not discussed here. 


2.2.3. The Iteration Count Coloring 


One of the most intuitive ways to specify the color of a pixel is to choose the 
color proportional to the number of iterations N. The Iteration Count coloring 
is illustrated in Figure 2.1. The coloring function simply returns the number of 


elements in the truncated orbit divided by 20. The coloring is applied to the 
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system 2.3 with c = 0.5 + 0.257 and p = 2. The bailout value was chosen to 
be M = 10 and the initial value z) was set to the value corresponding to the 
pixel position. The real and imaginary parts vary between -1.333 and 1.333 in 
the figure. The color palette function uses the RGB color palette shown in figure 


2.2. The color index function 
I(u) = k(u— uo), 


where k = 2.5 and uo = 0, was used. 

This example fractal is used throughout this work to illustrate different col- 
orings. The coefficient k is used to adjust the colors in the image to a suitable 
range. The offset uo is used to adjust the relative location of light and dark 
areas. The outside areas of fractals are of most interest for artistic purposes. 
This fractal consists only of outside points and is consequently well suited as an 
example. The ideas presented can easily be adapted to many other divergent 


fractals as well. 


2.3. Iteration Levels and Boundaries 
The following definition assumes that a function f and a bailout value M are 
given. 


Definition 2.5 (Iteration Levels and Boundaries). For n > 1, the nth iteration 


level is the set 
Ly = {zo € C: |f""(z0)| < M < |f"(z0)|}- 
The nth iteration boundary is the set 
By = {ee eCz|f" *(2)|= ™M). 

For n = 0, define 

Lo = { f (20) 72 € Ly} 
and 

Bo = { f (20) PZ0°E B,}. 


The iteration levels correspond to the areas colored with same color in Figure 
2.1. The boundaries lie between the iteration levels. A point in Lj) U By exceeds 
the bailout M after N iterations. 
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A, 


Figure 2.1: Illustration of the Iteration Count coloring. 


2.4 Desirable Properties of Colorings 


The following subsections present criteria to evaluate the goodness and usefulness 


of a coloring. 


2.4.1 Performance 


The fewer calculations needed to determine the color of a pixel, the better. At the 
moment of writing, rendering times for fractal images on a normal PC are usually 
manageable. For example, rendering the final image constructed in Appendix A 


using both cores of a 2.17MHz Intel Core 2 processor took 8 minutes and 28 
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0 


Figure 2.2: The RGB color palette used in Figure 2.1. The colors are indexed 
from 0 (left) to 1 (right). 


seconds. The image was rendered using resolution 2000 x 1225 and normal anti- 
aliasing settings. 

A fractal animation sequence can contain thousands of frames and each frame 
must be rendered separately. Consequently, rendering times for fractal animation 
are orders of magnitude greater than rendering times for images. Rendering a 
short animation sequence can take from hours to days and thus performance 


remains to be a problem. 


2.4.2. Smoothness 


It is often desirable that the colors of a fractal image change smoothly, that is, 
without noticeable edges. In addition to any discontinuities that are very easy 
to see, the human eye also identifies discontinuity in the rate of change of colors. 
Continuity is usually sufficient if the discontinuity in the derivative is small. 
However, smoother appearance requires the coloring function to be continuously 
differentiable. 


Definition 2.6 (Smooth Coloring). A coloring u is smooth in z if u(Or(zo)) € 
C1(C), that is, if both u(Or(z0)) and 


du(Or(z0)) 
dz 


are continuous [5, Remark 3.1.2, p. 251]. 


Colorings are discontinuous in the neighborhood of the periodic points P of 


a fractal set. Define 


Pr={e CC: Spe Pe |2=p< Tr}. 
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The set P, includes a neighborhood of radius r to each periodic point in P. In 
this work, a coloring is said to be differentiable if it is differentiable in the set 
{z €C: |z| < M}\P, for a small r > 0. Thus the neighborhood of the periodic 
points is ignored when considering the differentiability of a coloring. 

Discontinuity usually appears at the iteration boundaries because the number 
of iterations changes at the boundary. The smooth iteration count and interpo- 
lation can effectively be used to eliminate discontinuities, as will be discussed in 
Chapter 4. 


2.4.3 Degrees of Freedom 


A parameter is a variable in the coloring function that affects the visual appear- 
ance of the coloring. Parameters offer degrees of freedom for adjusting the visual 
appearance of the coloring. 

Parameters that give smooth coloring only for specific (e.g. integer) values 
appear frequently in applications. Real parameters that continuously affect the 


coloring are especially useful for creating motion in fractal animation. 


Chapter 3 


The Smooth Iteration Count 


Coloring 


The smooth iteration count is a smooth equivalent of the iteration count. Both 
attain integer values on the iteration boundaries, but the smooth iteration count 
varies smoothly between the iteration boundaries whereas the iteration count 
changes in a stepwise manner. 

The smooth iteration count is useful for creating smooth and soft color tran- 
sitions in fractal images. Its decimal part can also be used for interpolation in 


order to introduce continuity and smoothness to various other colorings. 


3.1 Simplification 
Consider the system 2.3. The approximation 
Gao ater; (3.1) 
is valid for |z?_,| >> |c|. This gives the formula 
Pe ee (3.2) 
Denoting zp = roe’ and applying 3.2 iteratively gives 
Zn = vb elPF0. (3.3) 


or 
l2n] = ro". (3.4) 


14 
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The bailout inequality 2.1 states that M < |zy|. The function |f(z)| = |z|? 
is strictly increasing for |z| > 1. Assuming M > 1 and applying the formula 
lzw| = |znw_1|? to the inequality |zy_i| < M gives |zy| < M”. Thus the last 
iterate |zy| fulfills the inequality 


M < |zn| < M?. (3.5) 


Figure 3.1 shows the iteration count coloring applied to the system 3.2. The 
figure shows the square bounded by the points —4—4i and 4+ 42 in C. The black 
disk |z| < 1 in the middle of the image consists of inside points. Consequently, 
only the distances rp > 1 are relevant for this system. 

For the system 2.3, it can be shown that all orbits O(z) escape for |z| > 
0.5 + \/0.25 + |e] [1, p. 258]. Assuming |c| < 2, all orbits for |z| > 2 escape so it 
is reasonable to confine rp to values greater than 2. 

The approximation 3.1 offers an accurate enough simplification of the sys- 
tem 2.3 for large ro. Most importantly, the nth iterate can be expressed in a 
very compact form according to equation 3.3, whereas the system 2.3 leads to 


polynomials of increasing complexity that are difficult to work with. 


3.1.1 Extending the Results 


The Continuous and Smooth Iteration Count colorings are deduced for the sim- 
plified system 3.2. It will be shown in Section 3.3.2 that the Smooth Iteration 
Count coloring is almost smooth also when applied to the system 2.3. 

When applied to the system 3.2, the Continuous and Smooth Iteration Count 
colorings are independent of the argument 6p of the initial point z9 = roe’. Thus 
it is sufficient to consider the continuity and smoothness with respect to the initial 
radius 1. 

Equation 3.3 allows the representation of a coloring u(ry) as a function of ro, 
namely 

uty): = u(r?” ). 
Expression u(re” ) is often used for deducing results for the system 3.2. However, 
when extending the results from the system 3.2 to the system 2.3, it is necessary 
to use u(ry) because the simplification 3.1 is likely to be valid only for the last 


iteration. For this reason, also the discussion of continuity and smoothness of 
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the colorings operate on u(ry). The requirement ro > 2 is also replaced by the 


requirement Try > 2. 


Figure 3.1: The iteration count coloring applied to the system 3.2. 


It is noteworthy that for any outside point zo, its iterates 21, 22,... start to 
grow at an increasingly fast rate as soon as they escape outside the area where 
system 2.3 exhibits chaotic behavior. Formally, for an arbitrary bound R there is 
a number n such that the iterates z,,Z,11,... are outside the disk |z| < R. The 
requirement 79 > 2 and the assumption made in the next chapter that ro is large 
are valid for the iterates z,, 2n41,... When extending the results to hold for the 
system 2.3, the iterates z9,...,Z,—1 can be omitted because the Continuous and 


Smooth Iteration Count colorings are functions of the magnitude rj of the last 
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iterate only. Thus the conclusions made are also valid for outside points zo for 
which 79 < 2. 


3.2 The Continuous Iteration Count 


A continuous but not Smooth Iteration Count coloring is first presented. This 
coloring exhibits noticeable edges at the iteration boundaries. It is presented in 


order to illustrate and emphasize the requirement for smoothness. 


Definition 3.1 (Continuous Iteration Count Coloring). The Continuous Itera- 
tion Count coloring is the function 


M? —-rn 


u(ry) =N+ (3.6) 


According to the bailout inequality 3.5, the magnitude ry of the last iterate 
varies between M and M?. The decimal part of the continuous iteration count 
is given by the quotient 

M?—-r N 
MP — M- 
This expression can be interpreted as linear interpolation in the interval |M, M?]. 


It attains the value 0 when ry = M? and 1 when ry = M. 


3.2.1 Continuity and Smoothness 


The following theorem describes the continuity and smoothness of the Continuous 


Iteration Count coloring when applied to the system 3.2. 


Theorem 3.1. The Continuous Iteration Count coloring 8.6 attains values equal 
to the iteration count at the iteration boundaries. It is continuous but not smooth 


with respect to ro. 


Proof. Continuity on the iteration levels follows from the continuity of the ex- 
pression 3.6. The following shows that it is continuous also in the neighborhood 
of the iteration boundaries. 


Consider a point 29 € By. The definition of the iteration boundary gives 
|zn-1| =M 


or, combining this with 3.4, 
lzw| = M?. (3c) 
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Substituting to 3.6 gives 


u(M?) = 


This also shows that the continuous iteration count coincides with the iteration 
count at the iteration boundaries. 

Now consider another point 2) € Ly_, such that |Zy_,| = M +e for a small 
€ > 0. This orbit escapes outside the bailout M at the (N — 1)th iteration. 
Denote ¢ | 0 ase —Oande> 0. Substituting ry_; = M +e to 3.6 gives 


‘ MP -(M+e 
u(fy-i1) = N-1 Nt ) 
— N-1+1 
— aes 


as € | 0. Thus 3.6 is continuous in the neighborhood of the iteration boundaries. 


Using Equation 3.4, u(ry) can be expressed as 


pn 
u(r?) =N+ —— ; 
The derivative can now be calculated as 
du op” rent 
dro Me—M- 
for zp. For the point 29, the derivative is 
du pr-1pps 4 
dry | MP—M 


~1(p-1) 


The derivatives differ by a factor of pre so 3.6 is not smooth in the neigh- 


borhood of the iteration boundaries. 


Figure 3.2 illustrates the continuous iteration count as a function of the mag- 
nitude of the initial point ro. Figure 3.3 shows the Continuous Iteration Count 
coloring applied to the example fractal described in Section 2.2.3. Although the 
iteration boundaries are softer than in Figure 2.1, they are still clearly distin- 


guishable despite the continuity of the coloring. 
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Continuous iteration count coloring 


4.5 T T T T T T T T 
— lteration count N 
4 z --+» Continuous iteration count u 4 


Figure 3.2: The continuous iteration count u(re”) as a function of rp. The curves 


are calculated for M = 5 and p= 2. 
3.3. The Smooth Iteration Count 


The following discussion is based on the simplification 3.4. After the smooth 
iteration count expression has been deduced for the simplified system, the result 
is shown to be quite accurate for the original system as well. 


Applying Equation 3.4 to the bailout inequality gives 
M <r? < MP. (3.8) 
By taking a logarithm twice over this inequality, the bounds for N can be solved; 


nM < p¥Inro <plhnM 


In M In M 
in < Ninp <Inp+In = 
nT Inro 
1 In M 1 In M 
Spee I See (3.9) 
Inp In7ro Inp Inro 


The number of iterations N is an integer and its value changes by one when 


ro crosses an iteration boundary. The function 


1 In M 
= 1+ — In —— 3.10 
v(ro) ' In p . Inro ( ) 
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Figure 3.3: The Continuous Iteration Count coloring applied to the example 


fractal. 


attains positive real values for 2 < rg <M. Because N = N(ro) is bounded by 
the inequality v(r9) — 1 < N(ro) < v(ro), it is apparent that v(ro) must follow 
N(ro). Also, v(ro) = N(ro) holds at the iteration boundaries. This is illustrated 
in Figure 3.4. 


Using 3.4, 79 can be expressed as a function of the last iteration ry. Substi- 
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Smooth bounds of the iteration count coloring 
4.5 T T T T T T T 


T 
—— Iteration count N 
4y -—- LowerboundL 7 
: - Upper bound U 


Figure 3.4: The iteration count N is bounded by the upper bound U = v(ro) 
and the lower bound L = vu(ro) — 1. The curves are plotted for M = 5 and p = 2. 


tuting ro = re to u(r) and denoting u(ry) = v(re,) gives 


1 In M 
ulin) = In = 
np iInry 
1 pN In M 
ae In 
Inp lnry 
1 In M 
= 14+—(NInp+In— ) 
Inp lnry 
1 In M 
= N+1+—hn— 


Inp "nr N- 
Definition 3.2 (Smooth Iteration Count). The Smooth Iteration Count coloring 


is the function 
1 In M 


=N+14 ; 
u(rw) Inp ”inty 


(3.11) 


3.3.1 Continuity and Smoothness 


Theorem 3.2. The Smooth Iteration Count coloring 8.11 attains values equal 
to the iteration count at the iteration boundaries. It 1s continuous and smooth 


with respect to ro. 
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Proof. In analogy with the proof of Theorem 3.1, ry = M? for an arbitrary point 
zo € By. Substituting to 3.11 gives 


1 In M 
M?) = N+1 
uM) Inp InP 
1 In M 
= N+1 | —] 
ne “inM np) 
1 
= Nota 12 int 
Inp 
= N. 
All outside points on any of the iteration levels Ly, N = 1,2,..., are mapped 


by f% to the last iteration level Lo. Similarly, all bounds By are mapped to the 
last bound Bo. For any fixed outside point, the number of iterations N is finite 
so the mapping f is well defined and smooth because f is smooth. 

Also u(r) is smooth for a fixed N. Smoothness on any outside point zo € Ly, 
N = 1,2,..., contained in any of the iteration levels then follows from the 
smoothness of f%. 

The number of iterations N changes at the iteration boundaries. Because 
points in Ly, and By are mapped smoothly to Lo and Bj, respectively, the 
continuity of wu and oe needs to be shown only between Lp and By. 

Consider the points zo € B, and 2) € Lo such that r9 = M and fp = M +e. 
It was already shown that u(ry) = N for ry = M?, so u(r1) = 1. Substituting 
fo = M+e to 3.11 gives 

i In M 
Inp a In(M + ¢) 


uMt+e) = 0414 


1 
— 1+ —lnl 
Inp 


a | 


as € | 0. Thus 3.11 is continuous in the neighborhood of iteration boundaries. 
Differentiating 3.11 with respect to ry gives 


du 1 1 In M 


drn = In p In M/Inry ry(Inry)? 
1 


rvnlurylnp 


(3.12) 


Using the chain rule, the derivative with respect to rg can be written 


du du dry 
dry  drn dro 
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For ro = M, the last iterate is ry =r, = M”. Thus 


dry 
— = = pir} 
dro P 
and 
d MpP-t 
Be ao oe es (3.13) 
dro MP? ln MP lnp 
1 
=> ——., 3.14 
Min M Inp ( ) 
For 79 = M +e, the number of iterations is N = 0. The derivatives now evaluate 
to 
dry _ 
dro i 
and 
du 1 
dro  (M+e)In(M+e)Inp 
1 
Min M Inp 
as eé | 0. That is, “ is continuous between B, and Lo. 


’ dro 


The Smooth Iteration Count coloring is illustrated in Figure 3.5. The itera- 
tion boundaries are now indistinguishable from the figure. This is a consequence 
of the fact that the derivative i is continuous at the iteration boundaries for the 
Smooth Iteration count coloring but discontinuous for the Continuous Iteration 
Count coloring. 

The N — 1 first iterations only increase the iteration count, whereas the 
decimal part of u(ry) only depends on ry. Thus the N — 1 first iterations 
are independent of the approximation 3.1, which is important only for the last 
iteration. The statement made above about the continuity of f’ holds for the 
function f(z) = 2”, that defines the simplified system 3.2, as well as for f(z) = 
z? +c. Theorem 3.2 can thus be extended to the system z,4; = 2? + c as long 


as the approximation 3.1 is valid for the last iteration. 


3.3.2. Error Estimation 


Theorem 3.3. Applying the Smooth Iteration Count coloring given by 3.11 to the 


system 2.2 results in a discontinuity between the iteration levels. the magnitude 
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Figure 3.5: Smooth Iteration Count coloring applied to the example fractal. 


of the discontinuity is bounded by 


1 In MP 


Dade a . 
Inp f In(M? — |c|) 


(3.15) 


Proof. Consider again the two points zo € By, ro = M, and 2, f97 = M +e. If 
the system 2.3 is iterated instead of the system 3.2, a discontinuity is introduced 
between the iteration levels. For zo, the magnitude of the first iterate is ry = 


|z6 +c| » M?, so N = 1. This gives the smooth iteration count 


1 In M 
Wr) = 2 n 
Inp In|zi+c 
1 In MP 
ee (3.16) 


Inp "in| +e 
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For 2), N = 0 and 


: 1 In M 
UNO) Se Inp ; In(M +¢) 
wae ee 2s x InM 
Inp mM 
= 1. (3.17) 


The magnitude A of the discontinuity is given by 


A |u(r1) — u(ro) (3.18) 


1 In M? 
= |14 In 7 —1| 
Inp = In|z +e 
1 In M? | 


n 
Inp In|zi+ce| 


The triangle inequality gives bounds for |z4 + c|, namely 
261 — lel S |zo + el S |26] + lel, 
or, substituting r9 = M, 
M? — |c| < |e6 tel < M? + cl. 


The magnitude of the discontinuity A is to be maximized on this interval in 
order to find an upper bound for it. Denoting x = |zj + c| and differentiating A 


gives 
dA(x) | 1 
dx  a«lmalp 


Assuming x > 1, the derivative is strictly decreasing and 3.18 attains its maxi- 
mum value at |z} + c| = M? — |c|. Thus the maximum error due to the approxi- 


mation 3.1 is 
1 In MP 


Rese 
me Inp " Tn( MP — cl) 


For the example fractal with parameters p = 2, M = 10 and |c| = 0.5590, 
3.15 gives the maximum error of Ajax = 0.00176. This is illustrated in Figure 
3.6. Error equal to 0 is colored white whereas error equal to 0.00176 is colored 
black. 

Suppose the palette function maps the unit interval [0,1] to, say, 500 discrete 


colors and the index function is J(u) = u. Then the colors are separated by steps 
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of 1/500 = 0.002 in u. This is greater than the magnitude of the discontinuity 
Amax = 0.00176, so the extension from the simplified system 3.2 to the chaotic 
system 2.3 indeed only leads to an insignificant error. In addition, the accuracy 


can be further improved by increasing the value of M. 


V 


Figure 3.6: The error in the smooth iteration count due to the approximation 
Sel; 


Chapter 4 
Branching Average Colorings 


This chapter introduces the Triangle Inequality Average and Curvature Average 
colorings. Their properties, similarities and differences are analyzed and, based 
on this analysis, a definition of branching average colorings is presented. Section 
4.5 introduces a new branching average coloring called the Stripe Average. Its 
use is exemplified in Appendix A. 

All colorings named above are examples of average colorings. Before contin- 
uing to the discussion of these colorings, the average colorings must be defined. 


Also, two different interpolation methods for average colorings are presented. 


4.1 Average Colorings 


Given a truncated orbit O7(zo), introduce the notation 2” = {zj-m, Zi—-m41,---, 2i}- 
The notation Z7"(zo) = Z}” is used equivalently when the point zo used to cal- 


culate the iterates is important. The constant m is the number of iterates in Z}”" 


preceding z;. Average colorings calculate an average of a function t: C™*t! WH R 


evaluated at all points Z/"(zo) contained in a truncated orbit O7(zo). For a point 


zo € L; U B;, denote this average sum S;. Formally, 


S(Or(%))=—— YD tr (4.1) 


n=m+1 
where t, = t(Z/"(z9)). The function t is called the addend function and it is 
assumed to be a continuous function of elements in Z}”. 
Each element t(Z?") of the average sum requires m iterations that precede 2;. 


When applying average colorings to a fractal, the iteration levels Lo,..., Lm are 


Ze 
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excluded from the coloring calculations because the required number of preceding 


iterations are not available. 


4.1.1 Using Smooth Iteration Count for Interpolation 


Let i > m+1. Consider a point z € B;, 2 = roe’, and another point 2) € L;_1 
such that 2 = (ro + e)e™, where ¢ is small and nonnegative. The corresponding 


average SUMS are 
a 


for zp and 
$.-1(Or(4)) = ——_ > 27") 


for 2. For all elements 2, € Or(Z9), Zn — 2m as € | 0. It follows from the 
continuity of t that also t(Z/"(Z,)) — t(Z?"(z)). The sum S;_; can thus be 
calculated for Or(zo) as € | 0. 

A discontinuity appears for the average sums between the points zo and Zo 
because, in general, S; # S;_,. This discontinuity can be eliminated by interpo- 
lating between S;_; and S; in L;UB;. The decimal part d of the smooth iteration 
count is 0 at B;_; and converges to 1 in the neighborhood of B; that intersects 
L;. Thus d can be used as the interpolation coefficient. 


Using linear interpolation, define for z € L; U B; a coloring function 
uz(Or(z)) = dS; + (1 — d)Si-1. (4.2) 


For zo € B; as above, d= 0 and uy = S;_1. For 2 € Lj_1, d > l and uz > Sj_-1 
ase | 0. That is, uy is continuous. However, linear interpolation normally results 
in discontinuity of the derivative at the control points $;. This is illustrated in 
Figure 4.1. 

Smoothness can be achieved by using some suitable type of splines instead of 
linear interpolation. A good choice are the Catmull-Rom splines [4, p. 324-325]. 
For z € L; U B;, define 


ug(Or(z)) = Ho(d)S; + Hy(d)S;_1 + Ho(d)S;-9 + Hy(d)S;_3. (4.3) 


The polynomials Hp — Hz are listed in Table 4.1. The polynomials are chosen so 


that at d = 0, 
dus 1 


ae —_ 5 (Si-1 — Si-3) 
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and. at-d =I, 
dug 1 


dd 2 
That is, the tangent of the graph of u at d = 0 (or, equivalently, at d = 1) is 


(S; — S;_2). 


parallel to the line segment connecting the previous and the next control points. 


Notation Polynomial 
Ho 5(—d? + d3) 
Ay 5(d + 4d? — 3d?) 
Hy 5(2 — 5d? + 3d?) 
A $(—d + 2d? — d) 


Table 4.1: The Catmull-Rom spline weighting polynomials. 


The coloring 4.3 interpolates between the points S$; at d = 0 and S;_, at 
d = 1, whereas the points $;_3 and 5; are needed for specifying the derivatives at 
d = 0 and d = 1, respectively. The coloring ug is smooth also in the neighborhood 
of B;. It is here called smooth interpolation due to its smoothness. Linear and 
smooth interpolation are illustrated in Figure 4.1. 

The above discussion of average colorings is summarized by the following 


definition. 


Definition 4.1 (Average Coloring). Given an addend function t: C™ + R and 


average sums specified by equation 4.1, an average coloring is given either by 4.2 
or 4.3. 


For some fractals, the lack of smoothness for linear interpolation may create 
artifacts similar to those in Figure 3.3. For fractals generated with various values 
of c in system 2.3, these artifacts are subtle and almost indistinguishable (cf. 
Figure 4.3). 

For some other fractals, however, the artifacts are clearly visible. Figure 4.2 
compares the linear and smooth interpolation methods applied to a detail of 
a fractal set that clearly shows the artifacts. The Triangle Inequality Average 
coloring was used. The fractal is calculated in analogy with the example fractal, 
but system 2.3 is replaced by the system 


za 
ame ee a 
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Linear and spline interpolation 


2 T T T T T T T T T 
— Linear 
“++ Catmull-Rom 
4h Control points 
0 | 
$-1} 
Ob 
~3b 
-4 1 1 1 1 | 1 1 1 | 


2 2.1 2.2 2.3 2.4 2.5 2.6 2.7 2.8 2.9 3 


Figure 4.1: Comparison of linear and Catmull-Rom spline interpolation methods 


on an arbitrarily chosen set of control points. 


where d = —0.2013 + 0.56387 and c = 1.4686 + 1.2657. Using smooth interpo- 
lation removes the artifacts. On the other hand, the structure of the coloring 
becomes blurrier because the average sums are interpolated over four iteration 


levels instead of two. 


4.2 The Triangle Inequality Average Coloring 


The Triangle Inequality Average coloring was originally developed by Kerry 
Mitchell. Figure 4.3 shows the Triangle Inequality Average coloring applied to 
the example fractal. 

The original implementation used linear interpolation. It was designed for 
coloring the Mandelbrot set whereas a variation adapted for Julia sets is presented 
here. In the Julia set variant, the constant c in the discussion below denotes the 
seed of the Julia set (0.5 + 0.257 for the example fractal). In the Mandelbrot 
variant, c corresponds to the position of the pixel to be colored. However, both 


variants are based on the same idea and the two variations merely reflect the 
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~~ “a 


Figure 4.2: Comparison of linear (left) and smooth (right) interpolation. 


differences in the calculation algorithm for Julia and Mandelbrot fractal images. 
Their visual appearance is very similar. 


Consider the system 2.3. The nth iteration can be written 
Ce ae ee: 
Applying the triangle inequality to the right side of the equation gives 


llZn—al — tell S lena + el S lena + lel, 


or 


ll@n—al = lell S len] S l2nal + lel. 
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Figure 4.3: The example fractal colored with the Triangle Inequality Average 


coloring using linear interpolation. 


Introduce the notation 


Mn = ||Zn—al — Iel| 
and 
Mn = |%n-1| + le| 
for the upper and lower bounds of |z,|. Note that my, = m,(Zn-1) and M, = 


M,,(Zn-1). These bounds can be used to define a quotient in the interval [0, 1]. 


Definition 4.2 (Triangle Inequality Average Coloring). Define the addend func- 


tiont: C2 R as 


= \Zn| — Mn 


(Zs) ie 


(4.4) 
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Using this addend function, the Triangle Inequality Average coloring is given by 
Definition 4.1. 


Consider the two cases when the term |z,,| reaches its upper and lower bounds. 


1. Ifarg2?_, = —arge, 2?_, and c point at opposite directions in the complex 


plane and the lower bound is reached. Formally, |z?_, +c] = ||z?_,| — |el|. 
2. If arg z?_, =arge, z>_, and c are parallel and |z?_, + c| = |z?_,| + |cl. 


Thus t(Z}) varies as a function of arg z?_, between 0 when arg z?_, = — argc 


and 1 when arg 27_, = argc. 


4.2.1 Behavior for Distant Starting Points 


The following discussion describes the properties of an addend t, under the 
assumption that the magnitude of the starting point z is large so that the 
simplified system 3.2 is accurate for the n—1 first iterations. The addend function 
t of the Triangle Inequality Average coloring is based on the effect of c on the 
argument of iterates. Thus the system 2.3 is used for the nth iteration whereas 
the system 3.2 is used for the n — 1 first iterations. 

Denote c = r,e%. Assuming r, < |zo|, the difference in the bounds M,, and 


Mr can be simplified to 


Mn -— Mn = lenal+Pe— |l@n-al — Pel 


= 2r. 


because |z?_,| > r-. The nth addend t, can be expressed as 


Tn — ™n 

Mn -— Mn 
2 2 
ry My 


2r.(T, + Mn) 


th = 


The square of r,, can be written 


, = al 
= (2h_y+c)(2h_, +0) 
= (Feige 4+ pie) (2 ge Pare) 
mR tHE rll 4 emt) 
= Pt r2 + 2r?_sr.cos(On_1p — ). (4.5) 
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In analogy with 3.3, the (n — 1)th iterate is z,_1 = aa It follows that r,_; = 
re and 0-1 = p” 19. Substituting to equation 4.5 gives 


re = oP" 4 op? 4 Oph r, cos(p"Oo — Oe). 
Using the identities above, the nth addend t,, becomes 
r2 — m2 
2re(Tn + Mn) 
roe" +r24 ore cos(p"99 — 8.) — (reP" _ ore, +r?) 
2re(Tn + Mn) 

ore’ r.[cos(p"Oy — 8) + 1] 

arate tr, —T.) 

cos(p"O — 8.) +1 


a Tr nm Th 4.6 
ae ele LS relrG 46) 


1 1 
aS cos(p"A9 — 8.) + 5 (4.7) 


th = 


In the last step the approximations |z?" + c|/r?” = 1 and r./r?” = 0 have been 
used. 

Expressions 4.6 and 4.7 encapsulate the characteristic properties of the Tri- 
angle Inequality Average coloring. Firstly, 4.6 shows that the addend t,, is only 
weakly dependent of the magnitude rg of the initial point. Secondly, t,, varies 
sinusoidally as a function of of the initial angle 4). Each iteration increases its 
frequency by a factor p. Section 4.4.1 explains how these properties lead to the 
branching appearance that is characteristic to the Triangle Inequality Average 
coloring. 

Figure 4.4 shows t(Z})...¢(Zjz) and their mean value as a function of r; and 
0,. Note that r; and @; correspond to ro and 6 in the above equations because 
evaluating t requires two iterates. The iterates z,...z4 are calculated for the 
example fractal. The curves are in accordance with the above observations based 


on expressions 4.6 and 4.7. 


4.2.2 Behavior in the Neighborhood of Origin 


Theorem 4.1. The Triangle Inequality Average coloring is discontinuous at 


points that are mapped to 0. 


Proof. Assume that z,, = ce’’, where ¢ is small, is an arbitrary iterate that has 


been mapped to the neighborhood of the origin. The corresponding upper and 
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TIA sum terms and average as a function of the initial magnitude 


T T T T T T 
iE 
0.8 
fa BS EEE ESSE SS et i me Vf OS inet ele ee er cy ete eae de 
a 0.6 
wk 
0.4 | 
Esse! 
0.2- ah 
F) MERESERERETED Pere IRERERETEREST EREREREREEES DS IEREREREESERED TERREEXESEEEED NEAERREAREEADD Do Doe pe 
2 21 2.2 2.3 2.4 2.5 2.6 2.7 2.8 2.9 3 
r 
1 
TIA sum terms and average as a function of the initial angle 
7! 
= 
WwW 


Figure 4.4: The dependence of the addends ft; ...t, on 7, and 62 in the Triangle 


Inequality Average coloring. System 2.3 was used in the calculations. 


lower bounds are 
Masi =eP+ Te 
and 
Mnyi = le? —7-| =Te — €?. 
The denominator in expression 4.4 is 
Mri — Map = 28°. 


Using equation 4.5, the addend t,,; can be written 


Prt1 — ™Mnt+1 


tn41 = 
Maa = M44 
2 2 
Pray — ™Mnsq1 


2€? (Tayi + Mn41) 
EP + r? + 2ePr.cos(pO — 0.) — (e7? — 2e?r, + 1?) 
2€”(Tn4i + Mn41) 
2e?r.[cos(p@ — 6.) + 1] 
2€?(Tr41 + Mn41) 
r.[cos(p@ — 6.) + 1] 
(Thi +e — €?) 
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The magnitude of the (n + 1)th iterate is 
race = lee” tel re 


as e — 0. Thus 


r-|cos(po — 6.) + 1 1 1 
tn41 | we ) | =} cos(pa ~ 6,) + 5 


This limit is dependent on the angle @ so t,41 is discontinuous at z, = 0. 


Figure 4.5 shows the discontinuities of the Triangle Inequality Average marked 


with circles. 


Figure 4.5: The Triangle Inequality Average coloring applied to the example 


fractal. The discontinuous points that are mapped to 0 are marked with circles. 
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4.3 The Curvature Average Coloring 


The Curvature Average coloring was originally developed by Damien M. Jones 
for Ultra Fractal in 1999. It is illustrated in Figure 4.6. As the figure and the 


following analysis reveal, it is both visually and mathematically very similar to 


the Triangle Inequality Average coloring. 


Figure 4.6: The example fractal colored with the Curvature Average coloring 


using linear interpolation. 


The coloring is based on the idea of approximating the curvature of a curve 


defined by discrete points. The curve is in this case defined by the points of a 
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truncated orbit Or(zo). The approximation used is 
t(Z7) = arg(— 


where arg:_,,,)% denotes the argument of z in the interval (—7,7]. This also 


defines the averaging function in the sum 4.1. 


Definition 4.3 (Curvature Average Coloring). Define the addend function t : 


C2? R as 


en — &n-1 


t(Z7) = al (— 7,7] (4.8) 


2n—-1 — &n-2 
Using this addend function, the Curvature Average coloring is given by Definition 


af 


Unlike the Triangle Inequality Average coloring, the Curvature Average is 


continuous at points that are mapped to 0. 


4.3.1 Behavior for Distant Starting Points 


Suppose 7p is large and the simplified system 3.2 can be used instead of 2.3. 


Substituting z, = z>” to 4.8 gives 


Zn — 2n—1 
t( Ze) -= ‘| aitee 5 | 
( n) "" £n=1 — 2n—2 
D prt 
29 — *0 
=e arg TT pr-l pn—2 
“9 — 
n-1 1 = 
— larg 2 ?™ 1) 
_ er n— n—2 (py 
BB 1) 
N— 2a 2-1 lp, 
eee 2 (p D(gP (p-1) a 
= —1,0 m—2 ( —1) 
i ee 
nr—1(, N—2Zl(n " 
Assuming 28" ?-) > 2?" ®-) 55.1 gives 


nm—1 fm 
1a ea ee | 


Oop"""(p — 1) 


= |- 7+ Oop"""(p — 1) — 2a |. (4.9) 


where |a| denotes the largest integer less than or equal to a. 
The visual appearance of the Curvature Average coloring is in accordance 


with equation 4.9. Figure 4.7 shows t(Z3)...t(Z?) and their mean value as a 
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function of rp and 6). The iterates are generated for the example fractal. The 
addends are clearly almost independent of r9. The approximative expression 
4.9 to t(Z?) defines a continuous and periodic function that corresponds to the 
absolute value of a triangle wave function (cf. Figure 4.7). Its frequency becomes 
p times larger at each iteration. Thus the Curvature Average coloring exhibits 


the similar properties as the Triangle Inequality Average that lead to branching. 


Curvature sum terms and average as a function of the initial magnitude 


Curvature sum terms and average as a function of the initial angle 
7 4g . T T 


vA 


Figure 4.7: The addends t(Z3) ...t(Z?) as functions of rz and 62 in the Curvature 


Average coloring. System 2.3 was used in the calculations. 


4.4 Similarities and Differences 


Although the Triangle Inequality Average and Curvature Average colorings are 
based on different approaches, they appear to have visually similar structure. 
This is due to the fact that the addends exhibit similar behavior for distant 
iterates, as can be noted by comparing equations 4.7 and 4.9. 

As an iterate of an orbit escapes outside the area where the system 2.3 behaves 


chaotically, the following iterates diverge at an increasing rate. The conclusions 
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made above for distant starting points are valid for these iterates. The corre- 
sponding addends appear to give the average sum its characteristic appearance. 

This claim is supported by Figure 4.8. It shows the Triangle Inequality Av- 
erage and Curvature Average colorings with the addends corresponding to the 
last 5 iterations excluded from the average sum. In contrast to the distant it- 
erates, the remaining addends do not exhibit regular behavior described by the 


simplified system 3.2. Thus the characteristic branching structure is not present. 


Figure 4.8: The Triangle Inequality Average (top) and Curvature Average col- 


orings (bottom) with the last 5 addends excluded from the average sums. 


The addends for these colorings are periodic functions and each iteration in- 


creases the frequency by a factor p. The periodic function is the cosine function 
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for the Triangle Inequality Average and the absolute value of a triangle wave 
function for Curvature Average. The difference in smoothness with respect to 
the initial angle 0) emerge from these functions — the cosine function is continu- 
ous and smooth whereas the absolute value of a triangle function is continuous 
but not smooth at its maximum and minimum points. The points where the 


Curvature Average is not smooth are clearly visible as edges in Figure 4.6. 


4.4.1 Emergence of Branches 


The following definition summarizes the properties that lead to the branching 
appearance characteristic to the Triangle Inequality Average and Curvature Av- 


erage colorings (cf. Figure 4.3 and Figure 4.6). 


Definition 4.4 (Branching Coloring). An average coloring function is said to be 
branching if its addends t, = t(Z?") fulfill the following properties when applied 
to the simplified system 2.3. 


dn 


1. Fork =n,n—-1,...,n—m, the derivatives oo 


of an addend ty with respect 


to the magnitudes r;, are small. 
2. Addends are periodic functions of the argument of iterates. 
3. Each iteration increases the period by an integer factor p. 


Property 1 implies that the addends vary more noticeably with respect to 
the argument of iterates than their magnitude, and lines perpendicular to the 
iteration boundaries appear at the iteration levels. Property 2 guarantees that 
the lines are evenly distributed around the iteration levels. Property 3 causes 
the number of lines to increase by a factor of p from L;_; to L;. Interpolation 
of the average sums between two iteration boundaries create an impression that 
the lines divide into p branches. 

For the example fractal, p = 2 so each iteration divides the branches into 
two. This period doubling is a consequence of Property 3 above. Figure 4.9 
illustrates that the number of lines is doubled between each iteration boundary 
when moving towards the fractal set. The figure shows the average sums of 
the Triangle Inequality Average coloring applied to the example fractal without 


interpolation. 
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Figure 4.9: The average sums of the Triangle Inequality Average coloring applied 


to the example fractal without interpolation. 


4.4.2 Continuity and Smoothness 


Except for points that are mapped to 0, the Triangle Inequality Average color- 
ing is continuous and smooth in the iteration levels. Ignoring the effect of the 
small discontinuity in the decimal part of the smooth iteration count, it is also 
continuous but not smooth at the iteration boundaries when using linear interpo- 
lation. Using smooth interpolation guarantees smoothness also at the iteration 
boundaries. 

The same conclusions also hold for the Curvature Average coloring except for 


the lines that emerge at the extreme points of the triangle waves. On these lines 
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the coloring is continuous but not smooth. 


4.5 The Stripe Average Coloring 


A branching coloring similar to the Triangle Inequality Average and Curvature 


Average can be achieved with a more explicit approach. 


Definition 4.5 (Stripe Average Coloring). Define the addend functiont: C'> R 


as 


1 1 
iz) = 5 sin(s arg 2) + a (4.10) 
The constant s is called the stripe density. Using this addend function, the Stripe 


Average coloring is given by Definition 4.1. 
Using the approximative identity 6, = p"6o, the addend function becomes 


1 1 
(eas 5 sin(sp”09) + 5 


For s = 1, the addend function is equivalent to 4.7 after a translation by 6. with 
respect to 09. Consequently, the Stripe Average coloring fulfills properties 1-3 
of Definition 4.4. Its behavior is equivalent to that of the Triangle Inequality 


Average for distant iterates. 


4.5.1 The Stripe Density Parameter 


As its name indicates, the stripe density parameter s can be used to adjust the 
density of stripes. Large values increase the variation frequency with respect to 
6) and thus result in denser stripes. Non-integer values result in discontinuity. 
Figure 4.10 shows the Stripe Average coloring using s = 3 (top) and s = 10 
(bottom) applied to the example fractal. Comparison with Figure 4.3 shows 
that the branching structure is similar but stripes are denser in Figure 4.10. 
Bailout M = 107° was used in both images. 

It is noteworthy that increasing the bailout value has a similar effect of creat- 
ing denser stripes. An increase to a power of p in the bailout adds one iteration 
and a corresponding addend to every average sum. This addend has a frequency 
that is p times larger than the frequency of the preceding addend. However, 
numerical accuracy sets a limit to how large bailout values can be used. Thus 
a higher stripe density can be achieved even for small bailouts by choosing the 


value of s appropriately. 


CHAPTER 4. BRANCHING AVERAGE COLORINGS 44 


Figure 4.10: The Stripe Average coloring, using the stripe factors s = 3 (top) 
and s = 10 (bottom), applied to the example fractal. 


4.5.2 Behavior in the Neighborhood of Origin 


For z, = ce’, the expression 4.10 becomes 
1 1 
t(Zn) = =sin(s0) + = 
(2%) 5 sin(s y+ 5 


This is independent of ¢ but dependent on 6. Thus the Stripe Average coloring 
is discontinuous at points that are mapped to 0. The discontinuities are clearly 
identifiable in Figure 4.10. 

Figure 4.11 shows the Stripe Average applied to the example fractal with the 


6 last addends excluded from the average sum. In analogy with the Triangle 
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Inequality Average and Curvature Average colorings, the branching structure 


disappears as the addends corresponding to the distant iterates are excluded. 


Figure 4.11: The Stripe Average coloring, using s = 3, applied to the example 


fractal. Six last addends were excluded from the average sums. 


Kapitel 5 
Svensk sammanfattning 


Malsaéttningen med detta arbete var att studera matematiska egenskaper hos 
nagra algoritmer som anvands for att farglagga fraktalbilder. 

Betrakta funktionen f(z) = z? +c som definierar ett dynamiskt system Zz = 
zP_, +c. Beroende pa vardena av zp och c kan iterationerna z; = f(z), Z2 = 
f(<1),... strava till oandligheten eller oskillera periodiskt mellan ett andligt antal 
punkter. 

Den klassiska algoritmen for berakning av divergenta fraktaler kan utvidgas 
pa foljande satt genom att infora index-, palett-, och farglaggningsfunktioner. FGér 
en punkt zo € C, berakna iterationerna z; = f(z0), z2 = f(z1),.--,2n = f(zn-1) 
tills |zv| > M. Konstanten M definierar mangden |z,| < M i C dar fraktalen 
beraknas. Om villkoret |z,,| > M inte uppfylls for nagon iteration 2, 22,..-, ZNmax 
sa betraktas punkten hora till fraktalmangden. Konstanten Nya. definierar max- 
imala antalet iterationer. Detta arbete koncentrerar sig pa omraden dar |z,,| > 
uppfylls. 

Farglaggningsfunktionen avbildar mangden av iterationerna {2Zp, 21,..., 2} 
pa ett reellt tal. Detta tal avbildas sedan av indexfunktionen och palettfunktionen 
till den slutliga RGB-fargen. Dessa operationer utf6rs for varje pixel. [terations- 
gradnser forekommer mellan punkter som kraver olika antal iterationer f6r att 
avbildas utanfor omradet |z,| <M. 

Férutom att krava sa lite berékningsoperationer som mdjligt ar det ofta 
onskvart att farglaggningen ar glatt, dvs. att den ar kontinuerligt deriverbar 
i C. Sadana parametrar som ingar i farglaggningsformeln och kan anvandas for 


att justera farglaggningens utseende ar ocksa nyttiga. 
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Farglaggningsfunktionen for den klassiska iterationsraknaren kan definieras 
som 


uU=kN 


dar N betecknar antalet iterationer som behévs fér att |zy| > M och k ar en 
positiv konstant. Farglaggningens varde andras diskontinuerligt vid iterations- 
granserna och halls konstant mellan dem. 

Den kontinuerliga iterationsrdknaren ar en kontinuerlig men inte glatt 
farglaggning och den ger klart synliga granserna dar iterationsraknarens varde 
andras. Den glatta iterationsradknaren ar glatt 6verallt. Dessa farglaggningar an- 
tar varden lika med iterationsraiknaren pa iterationsgrénserna men varierar kon- 
tinuerligt mellan dem. Iterationsgranserna kan inte urskiljas fran en fraktal som 
farglagts med den glatta iterationsraknaren. En jamforelse med den kontinuerliga 
iterationsrankaren visar varfor glatthet ar en viktig egenskap hos farglaggningar. 

Den glatta iterationsraknaren ar dock glatt endast for systemet f(z) = z? 
medan den visar diskontinuiteter pa iterationsgranserna for det kaotiska systemet 
f(z) = 2 +c, c 4 0. En Ovre grans for diskontinuitetens storlek kan beraknas 
med formeln 

Rs 1 o In MP . 
Inp In(M?P — |c|) 
I praktiken ar diskontinuiteten obetydligt liten for IZ = 10 och blir 4nnu mindre 


(5.1) 


for storre varden av M. 
Medeltalfarglaggningar ar en familj av farglaggningar som utnyttjar 


medelvardessummorna 


1 a 
Si 205. Zixe 2523) = ian oe ty: (5:2) 
n=m+1 


Addenderna t, beraiknas med en sa kallad addendfunktion t : C™*! »% R 


som ar en funktion av z, och eventuellt ett antal foregaende iterationer 


Zn—1; 2n—25+++,2%n—-m- Addendfunktionens egenskaper bestammer farglaggningens 
utseende. 

Medeltalfarglaggningarna beraéknar summorna S;_; och 5S; och interpoler- 
ar mellan dessa. Decimaldelen av den glatta iterationsraéknaren anvaénds som 
interpolationsvariabel. Linjér interpolation ger vanligen goda och nastan glat- 
ta resultat, annars kan man anvaénda Catmull-Roms kurvor for att fa en glatt 


farglaggning. 
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Triangelolikhet- och kurvatur-medelvardesfarglaggningarna ar exempel pa 
medelvardesfarglaggningar. Addendfunktionen for den sistnémnda ar 


an — &n-1 


(5.3) 


tans en—1; Rad) =. al S (— 7,7] 
2n—-1 — %n-2 


dar arg:_,,;% betecknar argumentet av z pa intervallet (—7,7]. For 


triangelolikhet-medelvardesfarglaggningen ar addendfunktionen 


2n| — Mn 
LF, Se) = Fen] = Min =i (5.4) 
dar my = ||2n-1| — lel] och My = |2n_1| + lel. 


Dessa farglaiggningar ar kontinuerliga. Om Catmull-Roms kurvor anvands 
for interpolation sa ar triangelolikhet-farglaggningen ocksa glatt. Kurvatur- 
farglageningens farglaggningsfunktion innehaller absolutbelopp och ar darfoér inte 
glatt. 

Dessa farglaggningar har en likadan forgrenande struktur. Denna struktur 
foljer ur foljande egenskaper hos addendfunktionen. 


1. Addendernas derivator a med avseende pa iterationernas storlek rz, ar 


sma. 
2. Addenderna ar periodiska funktioner av iterationernas argument. 


3. Varje iteration 6kar pa perioden med en heltalsfaktor p. 


En medelvardesfargliggning som uppfyller dessa egenskaper  kallas 
férgrenande medelvardesfarglaggning. Egenskap 1 medfér att farglaggningen 
varierar tydligare vinkelratt mot iterationsgranserna och bildar linjer mellan 
dem. Fran egenskap 2 féljer att linjerna ar jaémnt fordelade. Egenskap 3 medfoér 
att antalet linjer 6kar p-faldigt. Pa grund av interpoleringen ser linjerna ut att 
forgrenas till p grenar mellan varje iterationsgrans. 

For stora varden av z, ar triangelolikhet-medeltalfarglaggningens addend- 
funktion approximativt 


1 1 
5 sin(arg Z,) + = (5.5) 


5: 
Genom att infora den sa kallade randparametern s fas addendfunktionen 
1 1 

t( Sn) = si n) + =. 

(24) ; sin(s arg Zp) + 5 

Denna addendfunktion definierar rand-medelvdrdesfarglaggningen. Randparame- 


tern s bestémmer randtatheten hos denna farglaggning. 


Appendix A 


Using Colorings on a Fractal 


Image 


This Appendix illustrates the construction process of a fractal image. It also 
exemplifies some possible uses of the Smooth Iteration Count and Stripe Average 
colorings, and shows that the colorings discussed can be applied to other systems 
than the system 2,4, = 2? +c used in earlier examples. 

The image is constructed from multiple layers each containing an independent 
fractal image. Layers are combined together using merge modes that define how 
they interact visually [6]. The merge modes used in this example are described 
in Table A.1. 

A fractal set called Rudy’s Cubic Mandelbrot, developed by Mark Townsend, 


was chosen to be used for the image. It is defined by the system 
ne te zs +dz,+e 


where d is a constant. Both zp and c are set to correspond to the position of the 
pixel to be calculated. The bailout value 10'° and the exponent p = 3 was used 
for calculating the smooth iteration count. 

The image was created using Ultra Fractal version 4.03. At the moment of 
writing, the Smooth Iteration Count and Stripe Average colorings are available 
in the Ultra Fractal Formula Database in the file jh.ucl. The Ultra Fractal 


parameter file for the image is available in jh.upr. 
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A.1 Intermediate Images 


Figures A.1 - A.9 present a sequence of intermediate images and associated notes. 


ie ad 


Figure A.1: After trying various values, the constant d was set to —0.7198 + 
0.9111i. Interesting shapes shown on the right side of the figure were found in 


the location indicated on the left. 
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Figure A.2: Three layers using the Stripe Average coloring with stripe densities 
s = 2,8 = 4 and s = 8 were created. The opacities of the layers were set to 
100%, 50% and 33% so that 1/3 of the colors of the final result comes from each 
layer. This averaging technique has been proposed by Kerry Mitchell [7]. 

Stripe density s = 1 did not give enough structure whereas s = 16 resulted in 
too dense stripes. The figure shows the three layers averaged with s = 16 on the 
right and s = 8 on the left. 

The color densities and palettes were adjusted to give a pleasant variation of 
dark and light areas. This requires that the dark and light spots of the averaged 


layers partially coincide. 
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Figure A.4: Colors were added to the averaged layers. Analogous green and blue 
hues were chosen for each of the three layers in order to create an interesting 
variation of color. A new, partially transparent Stripe Average layer was added 


to highlight parts of the image with a color lighter than the background. 
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Figure A.5: The colors of the Stripe Average layers were further adjusted. The 
black color was replaced by transparency on the two topmost averaged layers. 
This made the final result lighter and removed gray artifacts from the background 
that emerged from averaging white and black. 

A layer using the Smooth Iteration Count coloring was added using the Hue 
merge mode. However, the layer was removed as it did not improve the appear- 


ance of the image. 
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Figure A.6: Two layers with a smooth disk-shaped color transition were added 
using the overlay merge mode. The first of the layers darkens the background 
tints. The second one, masked to only affect the highlight layer, transforms the 


highlights to create more variation in color. 
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Figure A.7: The size, position and color of the disk overlays were further adjusted. 
The size of the highlights was reduced. Minor adjustments in color were also 


made to other layers. 


Figure A.8: Another color scheme was made as well, with cold blue and green 
hues in the spiral. However, the author found the warm hues of Figure A.7 more 


appealing. 
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Figure A.9: Two layers using the Smooth Iteration Count were added to highlight 
the fractal structure. The first layer uses the hard light merge mode with opacity 
50% to add color to the background. The second layer uses the overlay merge 


mode to lighten the fractal structure highlights. 
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Merge mode 


Description 


Normal 


Directly returns colors from the layer. 


Multiply 


Multiplies the layer with the underlying layers. The 
result is always a darker color, thus darkening the un- 


derlying layers. 


Screen 


Multiplies the inverse of the layer with the inverse of the 
underlying layers, and inverts that again. The result is 
always a lighter color, thus brightening the underlying 


layers. Screen is the inverse of Multiply. 


Overlay 


Multiplies or screens the colors, depending on the color 
in the underlying layers. Creates color blending effects 


between the layer and the underlying layers. 


Hue 


Returns the hue of the layer, and the saturation and lu- 
minance of the underlying layers. Colors the underlying 


layers with the hue of the layer. 


Hard Light 


Multiplies or screens the colors, depending on the color 
in the layer. Emphasizes the dark and light regions in 
the layer, while the areas with medium brightness be- 


come transparent. 


Table A.1: Description of layer merge modes [8]. 
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A.2 Conclusion 


Tables A.2 - A.4 show the layer contents and their cumulative effect to the 
resulting image. Layer names, opacities and merge modes used to combine each 
layer to the underlying cumulative image are given above the corresponding layer 
images. The merge modes are described in Table A.1. 

The Stripe Average coloring is well suited for the layer averaging technique 
— the result has more interesting structure than a single layer colored using the 
Stripe Average. This technique also offers more freedom to adjust the final result. 

Whereas the Stripe Average coloring has its characteristic striped appearance, 
the soft color transitions created by the Smooth iteration Count coloring direct 
the attention of the viewer to the fractal structure. Combined with the use of 
transparency, it is useful for highlighting the fractal structure without affecting 


other areas of the image. 
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Layer Image Cumulative Image 


3. Stripe Average s = 8 (Normal, 33%) 


Table A.2: Layers 1-3 of the image. Layer name and number is shown above 
the corresponding images. Merge mode and opacity are given in brackets. The 


checkerboard pattern in the images indicates transparent areas. 
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Layer Image Cumulative Image 


6. Lighten background overlay disk 


( 


Overlay, 100% 


) 


5. Highlights mask (Mask for layer 6) 


(fF Sam 


Overlay, 100% 


Table A.3: Layers 4-6. 
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USING COLORINGS ON A FRACTAL IMAGE 


Layer Image 


Cumulative Image 


9. Highlight spiral overlay (Overlay, 100%) 


Te ee knots (Normal, 1007) 


oer a 


Table A.4: Layers 7-9. 
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